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Abstract 

Some exact results are given, that connect, for a general magnitude of 
the interaction term U, the repulsive and attractive Hubbard model, in 
terms of the eigenspectra and quantum diffusion properties. In particular, 
it is shown that, for some initial conditions, the quantum evolution cannot 
differenciate between the attractive and repulsive models. These results 
apply to both fermionic and bosonic models, in any dimension and for 
general filling, as far as the underlying structure is bipartite. 

1 Introduction 

Among the numerous models that have been studied for years to describe corre- 
lated systems, the most popular is probably the Hubbard model|IJ, parametrised 
by a first-neighbor hopping term (5 and on-site interaction energy U . The parti- 
cles, electrons or bosons, can hop between orbitals located at neighbouring sites 
of a d-dimensional lattice. In one dimension, the exact eigenspectrum, for spin 
^electrons, is known, thank to the well-known solution by Lieb and Wu 
In higher dimension, no exact solution has been found, despite the large efforts 
intended to fulfil this task||. Therefore, partial results are welcomed , as far as 
they can provide a better understanding of this model. Usually, one considers 
separately the bosonic and fermionic versions of this model, as well as positive 
(repulsive) or negative (attractive) U on-site interaction. 

In this paper, I derive exact results that connect, for a general magnitude 
of the interaction term U, the repulsive and attractive model, in terms of their 
eigenspectra and quantum diffusion properties. It applies to both fermionic and 
bosonic models, in any dimension and for general filling, as far as the underlying 
structure is bipartite. There is already a long history of works devoted to using 



all sorts of symmetries that simplify the problem, some of which being intrinsic 
(like the SU(2) spin symmetry in the electron hamiltonian) , others being related 
to geometric properties of the underlying lattice or graph (like the bipartiteness 
property), and/or to particular electron filling (electron-hole transformation at 
half filling) ||||. In section 2, a symmetry property of the spectrum is derived, 
whch is certainly already known (see the side remark in ref p|), although I 
could not find it previously published in the present form. To the best of my 
knowledge, the consequence on quantum diffusion (section 3) is original. In 
any case, the demonstrations are not technically difficult, but some results may 
appear, at first sight, rather counterintuitive. Let me first summarize these 
results in two steps 

(i) it will first be shown that on a bi-partite structure, for general filling and 
dimension, the spectrum for positive U is opposite to the spectrum for negative 
U. In other words, to any eigenenergy E for one model, it corresponds (in a 
one-to-one manner) the opposite energy (-E) for the other model. The relation 
between corresponding eigenstates is simple and will be explicited. 

(ii) It is then possible to compare the time evolution, with the same initial 
state, for the two models. It will be proved that, for some class of initial 
conditions, these quantum evolutions cannot be distinguished. 

It is this latter assertion which may appear at first sight as surprising. As the 
simplest example, consider the case of two electrons of opposite spins located 
on the same site of a one dimensional chain at t = 0. If one computes the wave 
packet at any positive time, one finds that the weight (probability) of any ket in 
the tensor product basis is the same for U < and U > 0. The physical reason 
for that should appear clearly below. 

2 A symmetry property for the eigenspectrum 

Let us first prove the first part about the eigenspectrum property. We shall first 
construct a representation space £ for the many-particle Hilbert space, and map 
the many body problem onto a single-particle one in £ . 

2.1 Warm-up: the two-electron case 

To get familiar with this representation, let us describe the two-electron case for 
a one-dimensional chain. With two electrons of opposite spin, we consider one- 
particule states \l,a), where I runs on the sites and a = ±1 refers to the spin , 
which form a basis for the one-particle Hilbert space, and construct 2-particle 
states in the tensor product basis B 

B= | \M) = | Z,+)<g> | m,-) =| Z,m), I, me zj 

The representation space £ can safely be taken as a 2 dimensional square 
lattice, with a site (I, m) associated with the above ket | l,m). Since each 
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particle can hop only toward first neighbours in the chain, the non vanishing 
hamiltonian elements are restricted to first neighbours connections in £, that is 
along the edges of the square lattice. The standard Hubbard hamiltonian 



,+ 

la 



Cma + f/^C+ Q+C+Q- = T (P) + V(U), 



(1) 



<lm> ,cr 
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maps onto a one-particle standard tight-binding on a square lattice, with (3 
hopping terms, and a line of U on-site potentials on the main diagonal, 



H = (3^2(\l,m) (I + l,m\ + \l,m) (I - l,m| + \l,m) (l,m + l\ + (2) 



Note that, taking advantage of the periodicity along the diagonal (and con- 
structing Bloch sums labelled by k) , we can map this problem onto a family of 
one-dimensional problems with one impurity at the origin, and a fc-dependant 
hopping term. The standard Slater-Koster sum rule for this impurity prob- 
lem precisely corresponds to the Bethe-Ansatz condition for this two-electron 
problem Q]. 

Back to the 2 dimensional £ space, we note the mirror symmetry along 
the diagonal, which allows to split the solutions into fully symmetrical and 
anti-symmetrical states, living repectively in £s and £a- £s is a half space, 
containing the diagonal, defined for instance as the points {(l,m),l > m}, a 
given such point representing the symmetrized ket Ql, m) + \m, I)) / \/2. £a is a 
also a half space, excluding the diagonal, defined as the points {(l,m),l > m} 
, a given point representing the antisymmetrized ket (\l,m) — \m, I)) /\/2. It 
is then immediately clear that the antisymmetrical sector is [/-independent, as 
it is well known, and corresponds to the (addition) spectrum (without double 
occupancy) found in the problem of two electrons of equal spins. Indeed, in that 
case, Pauli principle imposes to restrict to the space antisymmetrical sector £a- 

The main points of interest here is that the two-particle problems has been 
mapped onto a single particle problem in a larger space with a U dependent 
subspace, and that the bipartiteness property of the original structure (the one- 
dimensional chain) is shared by the larger space (either £, £a and £s)- These 
properties remain for the general A^-body problem, as we shall see now. Since 
spin is conserved in the Hubbard hamiltonian, the iV-particle problem can be 
separated into independent problems where the number of particles of each spin 
is fixed. If particles live in a d-dimensional space J- ', the large space £ is just the 
euclidean sum of iV (a priori equivalent) copies of T . The spectral problem can 
be first studied (as will be done below) in terms of a single particle problem in 
£ , selecting then among the eigensolutions those satisfying the symmetries asso- 
ciated with particle spins. An alternative way is to first reduce £ by symmetry, 
and then solve the reduced hamiltonian. 
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2.2 Spectral symmetry 



Let us now derive our first result on the spectrum. We consider N particles, of 
specified spin, in a hypercubic lattice Z d , and construct iV-particle states |M), 
forming altogether a Hilbert space basis B represented by a Nd dimensional 
hypercubic lattice £ : 



B= I \M)= (g) \r i ,a i )\, (4) 

[ i=l,N J 

More generally, if the one-particule states lives on a bi-partite graph A (not 
necessarily hypercubic), with edges connecting kets with non vanishing hamil- 
tonian matrix elements, so do the /V-particule states. Indeed, let us first define 
a sign function s(r), taking the value +1 or —I according to which sub-lattice 
of A,T* belongs. In the hypercubic lattice case with unit edge, s(r) can for 
instance be defined as 

s(r) = (-l)^=i,^< (5) 

where xi are the d coordinates of ~f*. The sign function s(M) in the tensor 
product space is then simply defined as the product of the one-particule sign 
functions s(r) entering the tensor product, 

s(M) = [] s(n) (6) 

i=l,N 

Now, looking for non-vanishing elements (M'\H \M), with M' ^ M, it is 
clear that they occur only when one particule has jumped to a neighbouring 
site in A. This implies s(M') = —s(M), which proves that the A^-particule 
states live on a bipartite structure. The interacting part translates into N(N — 
l)/2 hyperplanes of dimension (Nd — 1), with U term as an on-site potential, 
corresponding to kets in B such that two particles share the same site in A. If 
these particles have equal spin, this manifold, that we shall call a [7-hyperplane, 
will be eliminated by a proper antisymetrisation (which amount here to build 
new kets as the difference between kets that are symmetrically related by this 
hyperplane mirror). Note that iteration of this process vith the whole sets 
of hyperplanes for equal spins is nothing but constructing Slater determinant 
states. The /7-hyperplanes intersect along hyperplanes of one dimension less, 
corresponding to states such that three particles meet on a A site, and so on, 
down to a single "diagonal" {/-space of dimension d, corresponding to states 
such that all particles share the same site. All these [/-hyperplanes play a 
prominent role in the bosonic case, while only the one with highest dimension 
occurs for the electrons problem. 

Since any symmetrization process done in £ amounts to generate the proper 
combination of sites related by the [/-hyperplane mirrors, it is not difficult to 
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show that the £ bi-partitcness property is shared by any of the properly sym- 
metrized or anti-symmetrized sets corresponding to the particle statistics. Those 
willing to consider a specific type of particle statistics, or to use more standard 
Fock spaces and second quantised operators may prefer to redo the following 
computations in that framework, which will eventually turn quicker. But I pre- 
fer here to analyze the spectral properties in the (much) larger structure £ , 
to stress that the properties derived here ( on spectral symmetry and quantum 
diffusion) are statistics independent. 

Let A(M) be the set of kets, different from \M) , and connected to it by H, 

A(M) = {\M') eB, M'^ M, and (M' \H\ M) ^ 0} , (7) 

and let \^(E,+)) be an eigenket of H(+\U\), with eigenenergy E, whose 
decomposition reads 

| *(£?,+)>= E°m|M) (8) 

MEB 

By projecting the Schrodinger equation onto the set of bra (M\, one gets a 
set of secular equations 

Yl a E M> = {E-f{M,U))a E M , (9) 

M'eA(M) 

where f(M, U) — (M \H\ M) is related to the number of multiply occupied 
1-particle orbitals in \M). The only property which we now use is that / is an 
odd function of U (it is in fact a number time U). Let us now consider the ket 

|*')- J2 s ( M )«m|M) (10) 

MEB 

where s(M) is the above defined sign function. Since M and A(M) live on 
two different £ subparts (for any M), \ *S?') clearly satisfies the set of secular 
equations 

P E a M , = (-E-f(M,-U))a M (11) 

M'eA(M) 

which proves that is an eigenstate of H(— \U\), with — E as its eigen- 
value: |\&') = \^(—E, — )). This correspondance can be carried out for all eigen- 
values in the spectrum, and we can write (Sp denoting the set of eigenvalues in 
the spectrum) 

Sp(H,U) = -Sp(H,-U) (12) 
At this point, two remarks arc worth pointing out 

(i) If one consider only changing the hopping term sign, one easily finds (a 
well known result) that 

Sp(H,-0,U) = Sp(H,l3,U) (13) 
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with, to a given eigenvalue E, |\&) as an eigenket for H(f3,U) and |^') for 

(ii) The hamiltonian H is the sum of a kinetic part T((3) and an interacting 
part V(U), In the £ space, as said above, this V{U) part is an on-site interaction, 
like in impurity models, whose value depends on the particular tensor product 
kets. The kinetic part takes a constant value whenever two kets are connected 
by H. This reminds of the one-electron so-called Anderson hamiltonian, but 
with two major difference: 

- In the Anderson hamiltonian, the diagonal terms are randomly distributed, 
while here the U dependant terms are heavily correlated. 

- One usually studies this Anderson hamiltonian onto lattices of physocal 
interest (in dimension 1,2 or 3), why £ is of high dimension. 

It is nevertheless worth noticing that the above symmetry on the two spectra 
apply as well to a standard Anderson hamiltonian, whenever all the random 
diagonal terms have their sign changed simultaneously. At this point, one should 
recall some recent work done on the two-electron problem (TIP) with Hubbard 
hamiltonian, in the context of possible interaction induced derealization S . If 
the (random, diagonal) potential distribution law is symmetrical around zero, 
the properties of the model ( like the eigenstates localization length) should 
not change if all diagonal terms have their sign changed. In such a case , one 
can predict that the interaction induced derealization around a given energy 
E (with Hubbard term U) should be identical to what occurs for —E and — U 

3 Quantum Diffusion 

We now compare quantum diffusion (evolution under the time-dependant Schrodinger 
equation) for both sign of U (attractive or repulsive interaction), and with spec- 
ified initial conditions. We first study the simplest case, where \^f(t = 0)) = 
\M) € B. As usual, the time evolution is computed by expanding \M) in the 
complete set of eigenkets £± = {\ty(E, ±))}. Let W (inverse of the above intro- 
duced aff) be the unitary matrix transforming the tensor product basis B into 
the eigenket basis £+. Let us stress that H being a real hamiltonian (which 
for instance would not be the case in the presence of a magnetic field), it is 
always possible to take W as a real matrix, which is now assumed. Now, with 
D = Card(B), we can write 

D 

\^(t = 0)) = \M)=J2w l M \^i), and 
l=i 

D 

|*(t)) = W i M eM-mt/h) I*,) (14) 
1=1 

where I labels the set of eigenstates of H{\U\). The componant of \^{t)) on 



G 



the ket \M') G 23 simply reads 



D rp 

(M'\^(t))=Y / (w l M ') W l M eM-iEit/h) (15a) 



1=1 



Now, with the same initial condition, let us compute the corresponding ket 
at time t , |\I/ _ (£)), for the hamiltonian H(— \U\). With the above discussed 
relations between the eigenvalues and cigcnkets, one easily obtains 

D rp 

(M'|*-(t)> = s(M) s (M')^(V z M ') W t M ' exp(iEit/h). (16) 

l=i 

Note that the product s(M)s(M') takes only two values, +1 if M and M' 
belong to the same sublattice, and —1 otherwise. So, in any case, we find a 
simple correspondance between \$?(t)} and \^~(—t)), 

(M'|*(t))=±(M'|*-(-t)) (17) 

Comparing both evolution for positive time, it comes 

\(M'\V-(t))\ 2 = (j^ (W^'f W t M eMiEit/t^ x 

(fl( W i M 'y W i M eM-iEit/h^J (18) 

= \(M'\y(t))\ 2 (19) 

So quantum diffusion, at any time t, is identical (in terms of probabilities 
with respect to any ket in the basis 23), for repulsive and attractive interaction 
of the same magnitude (\U\). This (surprising?) result is thus proved for any 
basis ket \M) , taken as an initial condition. It can be understood, on physical 
grounds, by recalling first that quantum diffusion is done at constant energy. 
The initial ket is delocalized in energy and its time evolution is influenced by 
low and high energy eigenket contributions. So, states close to the ground state 
energy for attractive U correspond to more or less gathered particles, while 
states with higher energy, eventhough we are in an attractive U model, corre- 
spond to more separated particles (this recalls the difference between standard 
bonding and antibonding states for one-electron states). So high energy states 
for attractive U influences the quantum evolution in a similar manner as low 
energy states for repulsive U, which qualitatively explain the above result. 

We now consider a more generic initial conditions for the quantum diffusion. 
Let us split £ (and therefore B) into two parts £a and £b , according to the 
bi-partite decomposition, (corresponding to a basis partition Ba and 23b), and 
write 

l*(* = 0)>= ]T &m,a|M)+ Y, b M,B\M). (20) 

MEB A MgBb 
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It is then easy to write the components of \^(t)) in the basis B: 



(M' |tt(t)) = J2 (W^'Y eM-iEit/h) 



1=1 



MeBa MgBb 



M 



The components of |* (i)) read accordingly 

d T 

(M' |*-(t)> - s(M') Y (^ M ') exp(iEit/h) x 



«(i4) ^ & M , A ^ M + *(£?) ^ & MiB ^' 
MeBa MeB B 



A I 



(21) 

(22) 
(23) 



It is clear that now, \(M' (t))f = |(M'|*(t))| 2 only if the initial ket 
|\Jj(t = 0)) lives on one sub-structure only (either A or B), with all components 
6m hi phase. This restricts the choice of initial conditions for which this (U, —U) 
parallel behaviour occurs to a large but definitively non generic set. Note how- 
ever that one is often inclined, owing to the linearity of Schrodinger equation, 
to restrict the numerical studies of quantum evolution to basis kets as initial 
conditions, which precisely belong to this non generic category. 
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